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Qj \ Abstract. 

"^ I Using exact expressions for the Ising form factors, we give a new very simple 

proof that the spin-spin and disorder- disorder correlation functions are gov- 

^ ■ erned by the Painleve III non linear differential equation. We also show that 

the generating function of the correlation functions of the descendents of the 
spin and disorder operators is a A^-soliton, N ^ oo, r-function of the sinh- 
Gordon hierarchy. We discuss a relation of our approach to isomonodromy 
deformation problems, as well as further possible generalizations. 
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1 Introduction. 

In the scaling limit, the Ising model is described by a free Majorana fermion 
i/jIx) with mass m = {T — Tc). However, the order and disorder operators 
(t(x) and /i(x) are non local in terms of these free fermions. This shows 
up in the very striking result by T.T. Wu, B.M. McCoy, CA. Tracy and E. 
Barouch [jl| for the correlation functions of a{x) and ii{x). Let r be the radial 
distance, then : 

</i(r)/i(0)>\ _ /coshix(s)\ „_/ 1 P^„, „,^^^X^2 .,„x.2 



< a(r)a(0) > ; - V sinh |x(.) ) ^"^ ['l X ^" H(^) " ^^^^ ^\ ) ' (1) 

where s is the scaling variable s = ^. x is a solution of the radial sinh- 
Gordon equation : 

(Py 1 dv „ . , /^ X 
Setting ?7 = e~^, this is equivalent to the Painleve III equation: 

v" fv'V 1 fv'\ , 2 1 



We present a new (as far as we know) very simple proof of it, based on the 
analysis of form factors. 

More generally, to compute exact correlation functions in 2D integrable 
quantum field theories remains a challenge despite recent progresses [@]. 
These progresses are based on a fine analysis of the algebraic Bethe anstaz. 
On the other side, it has recently been argued that the two-dimensional 
integrable QFT can be solved by only using their infinite quantum group 
symmetries P, ^. The simplest examples of such symmetries are provided 
by the Yangian symmetries of massive current algebras or by the quantum 
s/(2) loop symmetry of the sine-Gordon model. These symmetries act more 
simply on the asymptotic states and on the momentum variables. Therefore, 
they lead more directly to algebraic methods for computing field form fac- 
tors. The next question which naturally arises in this approach consists in 
knowing whether it is possible to reconstruct and to characterize, e.g. by dif- 
ferential equations, the correlation functions from the knowledge of the form 
factors in the spirit of the boostrap program [^]. In the following, we show 
by examining the simplest possible model, namely the Ising model, that it is 
not hopeless to expect a positive answer to this question. 



2 Ising Form Factors. 

As is well known 0, the Ising model is a Z2-invariant lattice model with 
nearest neighbour interactions between Z2 spins. It undergoes a second order 
phase transition at some critical temperature Tq- In the scaling limit near 
its critical point it is equivalent to a free Majorana fermion with mass m ~ 
(T — Tc). It possesses an infinite number of local integrals of motion in 
involution which have odd spins. The asymptotic states are the free fermions 
and therefore the scattering is trivial : 

^ = -1 (2) 

Among the fields of primary interest are the spin field a, its "dual" field 
called the disorder field /x, and the fermion 7/^. We refer to [0 for their lattice 
definition. 

Form factors are matrix elements of field operators. They satisfy alge- 
braic relations, usually called form factor axioms [||, |^, which depend on the 
locality of the fields and on the sectors to which the fields belong. Due to its 
Z2 symmetry, the sectors of the Ising model are labelled by couples of indices 
[a; 6] with a, 6 = 0, 1. In each sector [a; 6], the first index a refers to a rep- 
resentation of the group Z2 whereas the second index h refers to an element 
in Z2. The fields $[a;6](a;, t) in sectors [a; 6] satisfy the following equal-time 
braiding relations []: 

$[ai;6i](a;i;t) <l>[a2;fe2](x2;t) = (-1)"'^^ <I>[a2;b2](a;2;t) $[a^;fej](xi;t) ; Xi > Xs 

The identity operator is in the sector [0;0], the spin field a in [1;0], the 
disorder operator /i in [0;1], and the fermion '0, which is defined by the 

^ More generally, in a statistical model invariant under a finite group G, one naturally 
defines spin fields (t'', which take values in some representation p of G, disorder fields 
/ig, where the label g are elements of G, and parafermions '0^, which are defined by the 
operator product expansion -0^ ~ IJ-gc''. An element h £ G acts on the parafermions by 
TpP — !■ p(h)ipP. The parafermions ^^ span the sectors [p',g] which form representations of 
the quantum double 'D{G) of functions on G. The braiding relations of the parafermions 
are given by the universal i?-matrix of 'D{G): 

irg'M) ® 1)(1 ® rgl{X2)) = (1 ® P2{9l)i/l^^^-.{X2)){Vg'M) ® 1) ; ^1 > ^2 

The Ising model is a particular case of that. 



operator product expansion \1/ ~ /icr, is in the sector [1; 1]. In particular, the 
asymptotic particles, which are the fermions, are in the sector [1; 1]. 

In a more algebraic way, the couples [a; b] label the irreducible representa- 
tions of the quantum double I* (^2) of functions on the group Z2. In V{Z2), 
the universal i?-matrix acts on two representations [ai;6i] and [02; 62] as 
R = (— 1)'^^''^. The braiding matrices are therefore given by the universal R- 
matrix. Moreover, for the representation [1; 1] corresponding to the fermions 
we find R = —1. We therefore have the relation (which is a complicated way 
to write a simple result): 



R 



universal 



As expected from the general arguments ^ jTO) . 

Let us denote by \9i,---,9n) the n-particle states with energies Ei = 
mcoshOi and momenta Pi = msinh^j. The form factors for a field $[a;6] in 
a sector [a; b] are defined by : 

^S](^l,---,^n) = (0|$M](0)l^l,---,^n) (3) 

By crossing symmetry, all the matrix elements of $[a;6] are determined by 

^[a;b] ■ 

By Z2-symmetry, the form factors F}-^1 are non-vanishing only for a = 
n [mod 2]. 

The form factors axioms are usually written for local fields in the identity 
sector p, p. Their generalizations for non-local fields in the other sectors 
is easy to find. In ref [^ |T^, the generalized axioms were derived using 
quantum group symmetries. In the case of the Ising model, since all the 
sector are described by the Z2 x Z2 indices [a; b] and since the S'-matrix is 
S = —1, they reduce to : 

F[a;b]i^2,dl,d3,---,dn) = --^[a;6] (6*1 , ^2, ^'3, ■ ■ ■ , 6'„) 

Fgj(^i - z2n, ^2, ■ ■ ■ , ^n) = (-l)(i+'^)(i+'') Fg,(e,, ^2, ■ ■ ■ , ^n|4) 



<n) 



Res0^=e2-i7T F^ J 61,62,03,- ■■ ,9n) = (-1)" -1 F,,-, (6'3,---,6'„ 



\a+b 



.(n-2). 



[a;b] 



Here, we have used that a = n [mod 2] . Because there is no bound state, the 
form factors have no other pole in the physical strip, < Im^ < vr. 

In the sectors [a; b] with a + b = [mod 2], the form factors have no pole 
at all in the physical strip. They correspond to those of fields in the sectors 



of the identity and of the fermion. In the same way as in |TT[], the set of 
solutions to the equations (|) with a + b = [mod 2] are found to be in 
correspondence with the fields in the Neveu-Schwarz sector of the conformal 
Ising theory. 

The sectors with a+6 = 1 [mjod 2] are those of the spin and disorder fields. 
There, the form factors have a simple pole in the physical strip. This pole 
corresponds to the particle-antiparticle scattering. The minimal solution to 
eqs. (^) which possesses no zeroes in the physical strip is given by |^, [T^, |TT| : 



i^it(^i,---,^n) = n.tanhj^^). (5) 

These minimal form factors with n odd are identified as those of the spin 
field 0", and with n even as those of the disorder field /i. Eqs. (^ determine 
the minimal solution up to a multiplicative constant; we choose this constant 
to be one by a normalization convention; in particular (yu) = 1. 

As proved in |^ , all the solutions F^-^L to eqs. (^) can be written as prod- 
uct of the minimal solution times an auxiliary function, denoted P^"\9i, ■ ■ ■ , 6'„), 
which has no pole in the physical strip : 

Moreover, from eqs. (^), it follows that the functions P^"-* satisfy : 

P^^\e^-^27r,e2,■■■,e„) = P(")(^i,^2,---,^n) (6) 

P(")(^2+^Vr,02,^3,---,^„) = P("-'n^3,---,^„) 

In other words, P*^"^ are symmetric and periodic functions oi 9i,- ■ ■ ,9n and 
subject to the last constraints in eq. (^. It is easy to check that an infinite 
set of solutions to eqs. (@) is given by : 



<,...,.„}(%).=l,-,n = Qi?(%) •■■«(%) 



where si,---,sm, which label the solutions, are positive or negative odd 
integers. The functions Qs are : 



e ' 






As shown in ITT] , these solutions form a complete set of independent solutions 
to eqs. (^, and are in one-to-one correspondence with the fields in the 
Ramond sector of the Ising conformal theory. As a consequence, the form 
factors, 

j-,(2n+l) //) N _ p(2n+l) f a \ rpi'^n+^) ( a \ 

-^[l;0],{si,-,SM}^^i''j=l'-'2"+l - ^{s^,-,SM}y^i)^rnin K^j) 

-^[0;l],{ti,-,tQ}(^i)i=l,-,2n. = -P{ti,...,tQ}(%)-^mm(^i) (7) 

are identified with the form factors of the descendents as^^...,SM ^"^^ f^ti,--,tQ 
of the spin field a and disorder field fi. 

3 Ising correlation functions. 

We will only consider the two-point correlation functions. If Fi {6i, ■ ■ ■ , 9n) 
and F2 (Oi,- ■ ■ ,9n) are the form factors of two scalar fields $i(a;, t) and 
$2(2;, t), the Euclidean correlation function is then : 

($1(0:, t) $2(0,0)) = ^ -[ ^I'^^-mrcos^eA Ft\e„...,9^)Ft\e„---,9^, 

where r is the radial distance, r^ = x^ + 1^. 

Let us first consider the spin-spin and disorder-disorder correlation func- 
tions. Define : 

C±(r) = (/x(r)/i(0)) ± (a(r)a(O)) 

From the formula for the form factors, we have : 

Cdr) - I <f /:: n (*.—.) n -nh'^ (^) (9) 



We will study the sum of this series in the following section. 
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Consider now the correlation function of the descendents of the spin and 
disorder fields. Set : 

Since the operators a^'^^^ and fx^^''^ have non-trivial Lorentz spins, the corre- 
lation functions do not only depend on the radial distance but on x and t. 
We introduce the Euclidean coordinates z± = x ± it. From the expression of 
the form factors, we have : 



oo 



cty'^'^\x,t) =E ^ P n §e-f (-^-^--^^) ) Y[QsM)Q^M)m--^i'-^ 

n=0 ^' -^-^ j=l V ^^ / P,q t<j \ ^ , 

Here {sp} refers to multi-indices {sp^, ■ ■ ■ , Sp^^} and similarly {tq} = {t^j, ■ ■ ■ , tg^}. 
Let us introduce a generating function C^(zs) for these correlation functions. 
It depends on an infinite number of variables Zg with s = ±1, ±3, ■ ■ ■; we de- 
note z±i by z±. It is specified by : 



cF*.'(M)^n5^5^c,»w 



2s=0; |s|>3 



From the exact expressions (0) of the form factors of the descendents, we 
deduce : 



en.) - 1 ^/:; n i'^mi.)] n^^ (^) m 



where the "potential" X{9\zs) are given by : 



^ s=±l,±3,- 

In the next section, we show that the generating function C'^{zs) can also be 
expressed in terms of a solution of the sinh-Gordon equation. More precisely 
we will show that it is a r-function for the afiine sinh-Gordon hierarchy. 



4 Proof of the differential equations. 

The differential equation for the Ising correlation functions were originally 
proved by looking at the scaling limit of the lattice Ising model |]I| . See also 
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1^ . Another approach was developped by the Kyoto school which is based on 
the study of isonionodromy deformation problems ||T4]. We present another 
proof based on form factors. Our proof relies on comparing the generating 



functions Ca 



with a N —* oo limit of the iV-soliton r-functions of the 



afiine sinh-Gordon model. Therefore, we need to do a small detour into the 



afiine sinh-Gordon model. See ref. [^ for a recent study of this model. It is 
a Toda model over the afiine Lie algebra s/(2). It involves two fields, which 
we denote by (p and ^, whose equations of motion are : 



d^d^ (f) 
dudy ^ 



8M2sinh(20) 
SM^ cosh(20) 



Here M in the classical sinh-Gordon mass. The vacuum solution to these 
equations is (p^ac = and ^^ac = 2M^r^ with r the radial distance. Sub- 
stracting the vacuum solution, we have : 



dudu I 



8M^ sinh(20) 



8M^ 



cosh(20) - 1 



(11) 



where ^ = ^ — ^vac- Once the field is known, ^ is computed from it, if the 
boundary conditions have been specified. 

The afiine sinh-Gordon model possesses two r-functions t± which are 
defined by : 



r_ 



exp (-^) 
exp (-0) 



(12) 



The r-functions t± for the A^-soliton solutions are given by (See [|r^ and 
references therein) : 




p=l ki<k2<---<kp 



with 



X, = a,exp 2M^ 



Zs^J^^ 



Here Oj and /ij are the parameters of the N-sohton solutions. The variables 
Zs correspond to all the commuting flows of the affine sinh-Gordon hierarchy. 
The r-functions r^ are solutions of the quadratic differential equations : 

r±9,^9,_r± - 9,^r±9,_r± = -MV^ (14) 

These equations, which are the bilinear Hirota form of the affine sinh-Gordon 
model, are equivalent to the eqs. (|lT]) for the fields and ^. 

The claim is that C'^{zs) is a A^ — ;> oo limit of the functions r^ . The 
idea is to approximate the integrals appearing in the definition of C^ by 
Riemann sums. Therefore, let us set, 

m 
M = — ; Hi = -exp 

with A an arbitrary constant. Taking the limit A^ 




Finally, taking the limit A — > oo we exactly get the formula (0). Since the 
r-functions are solutions of the Hirota equations (|T3|) for any A^, a, and /Zj, 
it follows that C^ is also a solution. Setting : 

C^(^.)=expi(-eT0) 



then and ^ satisfy eqs. ([IT]). Integrating these equations and taking into 
account the fact that (p and ^ only depend on the radial distance proves the 
result of McCoy et al with the identification x = ~'P- The scaling variable s 
is defined as s = ^. 

Notice that, as it becomes to be familiar in quantum integrable models 
[0], the quantum correlation functions are the r-functions for a hierarchy of 
classical integrable non-linear differential equations. 

5 Miscellaneous remarks. 

• Relation with isomonodromy deformation problems. 
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In [Q , the Ising correlation functions were computed by using the Dirac 
equation for the fermions in an isonionodromy deformation problem. The fact 
that isomonodromy deformations are relevent to the computation of corre- 
lation functions can be understood from simple physical arguments. Indeed, 
monodromies of correlation functions reflect, and are actually equivalent to, 
the equal-time braiding relations among the fields. But the braiding rela- 
tions are scale invariant and therefore renormalization group invariant. In 
other words, the renormalization group transformations are isomonodromy 
deformations. 

Where do appear the isomonodromy deformations in our approach? We 
started from the algebraic equations defining the form factors, we solved them 
and then computed the correlations using the exact expression of the form 
factors. The point, making the relation with the isomonodromy deforma- 
tions apparent, is that the form factor axioms depend both on the matrices 
encoding the braiding relations and on the S'-matrix which also codes the 
braiding relations since it is given by the universal i?-matrix. Moreover, the 
mass scale does not show up explicitly in the form factor axioms. Therefore, 
computing the form factors is just an algebraic way to solve an iso-br aiding 
deformation problem. 

• Generalizations. 

Let us describe a way to generalize our approach to other situations. The 
proof that the correlation functions satisfy the afiine sinh-Gordon equations 
relied on their identification with soliton r-functions of this hierarchy. This 
identification is simply the remark that the r-functions and the form factors 
are built from the same kernel K{6i,6j) = tanh(^j — Oj/2). In the context 
of r-functions, this kernel is the expectation value of vertex operators from 
which the hierarchy can be reconstructed. 

This suggests the following procedure. First, assuming that the form 
factors have been determined (e.g. using algebraic methods inherited from 
quantum groups), find vertex operator for which they are the expectation 
values. 

Here the expectation values are in an auxiliary Hilbert space on which the 
vertex operators are acting. Then, reconstruct the hierarchy from the knowl- 
edge of the vertex operators (e.g. using Casimirs in the vertex operator 
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algebra). The correlation functions will then be soliton-type r-functions of 
the corresponding hierarchy. The simplest models in this approach are those 
with diagonal S'-matrices as we will describe elsewhere. 

6 Appendix: Fredholm determinant. 

Here, we show that the correlation functions can be written as Fredholm 
determinants. Indeed, using the identity. 




det 



\l^i + t^j J ^,J=l,..,n 



the r-functions C^{zs) becomes : 



oo /_i_Nn n4_p^ ™ 



c?w-i:^X:nU^A-(..i..) det 



n=0 "-• --^ i=l V" / YCOsh(^ 

By definition, this is a Fredholm determinant, 

C^{zs) = Det{l±K) 
for an integral operator with kernel. 




/,(»,»') = ,'-'^<''i=-> ' •'^■C'l^- 



V 2n cosh(^) V 2n 

It is interesting to note that all the dependence in the parameters Zg, as well 
as the mass scale, is concentrated in the potential X{6\zs)] i.e. the heart 
of the kernel, cosh"^((^ — 9')/ 2), is independent of all these parameters. 
The property, that all the coupling constants and all the renormalization 
group dependence manifest themselves only in the potential part of the ker- 
nel, appears frequently when dealing with Fredholm determinant in physical 
context, see e.g. 0]. 

Finally, notice also that since the square of the Ising model is equivalent to 
the sine-Gordon model at the free field point, this implies that the correlation 
functions of the sine-Gordon theory at that particular point can be written 
as the square of Fredholm determinants. 
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